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ABSTRACT 


The  natural  convection  heat  transfer  from  an  array  of  heaters  flush  mounted  on  a 
vertical  test  surface  in  response  to  both  step  and  periodic  input  power  has  been  investi¬ 
gated.  Baseline  surface  temperature  measurements  in  water  were  obtained  for  several 
steady  power  inputs,  ranging  from  0.2  \V  to  1.0  W  for  the  vertical  test  surface.  Also,  the 
effects  of  the  presence  of  a  parallel  shroud  near  the  test  surface  on  measured  temperature 
values  were  studied.  The  effects  of  periodic  variations  in  input  power  on  the  heater  ele¬ 
ment  temperatures  were  then  examined  at  selected  locations.  The  resulting  transient 
heater  temperatures  were  measured  for  various  amplitudes,  frequencies  and  average  lev¬ 
els  of  the  periodic  input  power.  These  transient  measurements  are  compared  with  a 
two-dimensional  numerical  simulation  of  the  relevant  fluid  flow  and  heat  transfer  proc¬ 
esses.  The  computed  heater  surface  temperatures  supported  the  trends  of  the  measure¬ 
ments  and  w'ere  within  the  expected  range. 
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I.  INTRODUCTION 


A  considerable  increase  in  the  volumetric  power  generation  rates  in  electronic  com¬ 
ponents  has  resulted  due  to  both  continued  growth  in  functional  capability  and  rapid 
decrease  in  feature  size.  The  future  enhancements  in  speed  and  size  reduction  of  elec¬ 
tronic  systems  are  constrained  by  available  heat  removal  techniques. 

The  performance  and  reliability  of  microelectronic  components  are  very  sensitive  to 
junction  temperature.  The  average  failure  rates  of  chips  increase  exponentially  as  the 
junction  temperature  increases  Bar-Cohen  and  Kraus  [Ref.  1].  In  recent  years,  due  to 
the  relatively  high  heat  transfer  rates  and  other  potential  advantages  such  as  high  reli¬ 
ability,  low  cost,  reduced  noise  and  simpler  design,  complete  immersion  of  microelec¬ 
tronic  circuits  in  dielectric  liquids  has  been  exanained. 

Relatively  few  studies  on  natural  convection  liquid  immersion  cooling  are  available 
in  literature.  The  heat  transfer  characteristics  of  flush  mounted  heat  sources  on  vertical 
surface  in  several  liquids  were  studied  by  Baker  et  al  [Ref  2).  Due  to  the  three  dimen¬ 
sional  effects,  large  increases  in  heat  transfer  were  noticed.  Both  flush  and  protruding 
arrangements  were  investigated  in  water  and  R-113  by  Park  and  Bergles  [Ref  3].  Cor¬ 
relations  for  natural  convection  in  a  vertical  channel  dissipating  uniform  heat  flux  from 
sides  were  presented  by  Bar-Cohen  and  Schweitzer  [Ref  4].  Flow  visualization  and  heat 
transfer  measurements  for  natural  convection  from  a  long  heated  protrusion  on  a  verti¬ 
cal  insulated  wall  of  a  rectangular  chamber  filled  with  water  were  carried  out  by  Kelleher 
and  Knock  [Ref  5).  These  measurements  were  supported  by  a  computational  investi¬ 
gation  by  Lee  et  al  [Ref  6). 
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Flow  visualizations  and  component  surface  temperature  measurements  for  natural 
convection  cooling  of  a  3  by  3  array  of  discrete  heated  protrusions  on  the  vertical  wall 
of  a  rectangular  enclosure  filled  with  a  dielectric  fluid  were  performed  by  Joshi  et  al 
[Ref  7],  In  other  experimental  studies,  the  heat  transfer  and  flow  characteristics  of 
protruding  heat  sources  on  a  vertical  surface  [Ref  8]  and  vertical  channel  [Ref  9]  were 
studied  by  Joshi  et  al.  Flow  visualizations  and  component  surface  temperature  meas¬ 
urements  were  carried  out  for  several  power  dissipation  levels  and  channel  spacings.  For 
the  smallest  channel  spacing,  a  significant  increase  in  the  component  surface  temper¬ 
ature  due  to  a  reduction  in  fluid  velocity  was  observed. 

Natural  convection  adjacent  to  a  three  column  array  of  fifteen  per  column  uniformly 
heated  components  flush  mounted  on  a  vertical  flat  plate  in  water  was  examined  by 
Gaiser  [Ref  10].  Component  temperature  measurements  both  without  and  with  a 
shroud  parallel  to  the  heated  surface  were  reported  and  experimental  correlations  were 
developed  for  different  power  levels  and  shroud  spacings.  A  significant  increase  in  tem¬ 
perature  due  to  reduction  in  fluid  velocity  for  the  smallest  charmel  spacing  was  observed. 

A  numerical  investigation  of  natural  convection  flow  and  heat  transfer  arising  from 
a  protruding  heat  source  on  a  vertical  plate  was  reported  by  Sathe  and  Joshi  [Ref  11). 
Numerical  predictions  were  obtained  for  a  wide  range  of  appropriate  Ra  and  Pr  numbers 
and  substrate  to  fluid  thermal  conductivity  ratios.  Computed  protrusion  surface  tem¬ 
peratures  compared  favorably  with  available  experimental  results  from  [Ref  8).  A  nu¬ 
merical  investigation  of  two  dimensional  natural  convection  flow  and  heat  transfer  from 
a  substrate-mounted  flush  heat  source  immersed  in  a  liquid-filled  square  enclosure  was 
conducted  by  Haukenes  [Ref  12).  Numerical  predictions  were  in  favorable  agreement 
with  the  experimental  data. 
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While  the  transport  characteristics  during  steady  state  have  been  investigated  for 
both  discrete  flush  and  protruding  heat  sources,  very  little  infonnation  is  available  dur¬ 
ing  transient  phases  of  operation.  In  many  applications  power  cycling  may  take  place. 
The  present  investigation  examines  both  numerically  as  well  as  experimentally  the  tran¬ 
sient  response  of  discrete  heat  sources  flush  mounted  on  a  substrate.  The  experiments 
were  conducted  with  the  same  test  surface  as  of  Gaiser  [Ref  10],  and  comparison  be¬ 
tween  numerical  model  and  experiments  were  made  for  the  heat  source  closest  to  the 
leading  edge  of  the  surface. 

The  experimental  studies  reported  here  examine  changes  in  the  natural  convection 
flow  and  heat  transfer  characteristics  of  a  three  column  array  of  45  flush  mounted 
heaters  on  a  vertical  test  surface.  The  power  to  the  bottom  heater  of  the  center  column 
oscillated  with  time  while  power  to  all  the  other  heaters  was  kept  constant.  Component 
surface  temperatures  were  measured  over  a  period  of  time  for  a  range  of  average  power 
inputs  from  0.2  W  -  1.0  W.  Various  amplitudes  and  frequencies  of  the  periodic  variation 
were  examined.  Also  reported  here  are  the  effects  of  the  presence  of  a  shroud.  Numer¬ 
ical  computations  were  performed  to  study  heat  transfer  from  a  single  heat  source. 

Specific  goals  of  this  study  were: 

•  To  measure  oscillations  in  heater  surface  temperatures  in  response  to  input  power 
fluctuations  for  a  range  of  average  component  power  levels,  amplitudes  and  fre¬ 
quencies  of  power  pulsations. 

•  To  examine  the  same  effects  in  the  presence  of  a  parallel  shroud  placed  next  to  the 
test  surface. 

•  To  examine  the  effects  of  power  pulsations  in  the  bottom  component  on  the 
steadily  powered  downstream  components. 
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•  To  numerically  investigate  the  two-dimensional  transient  natural  convection  flow 
and  heat  transfer  for  a  single  flush  heat  source  and  compare  the  results  with  the 
experimental  measurements  of  component  temperatures. 

This  study  advances  the  work  reported  by  Gaiser  [Ref  10]  and  Haukenes  [Ref  12] 
and  is  part  of  the  study  of  the  natural  convection  liquid  immersion  cooling  of  electronic 
equipment. 


4 


II.  EXPERIMENTAL  SETUP  AND  PROCEDURE 


Detailed  descriptions  of  the  original  experimental  equipment  design  and  con¬ 
struction  can  be  found  in  Gaiser  [Ref  10]  and  further  additions  and  changes  done  in  the 
experimental  setup  in  Haukenes  [Ref  12).  Only  a  brief  review  is  presented  and  new 
noteworthy  additions  and  changes  are  noted  here. 

As  seen  in  Figure  1,  the  assembly  consists  of  a  vertical  test  surface  with  a  total  of 
45  flush  mounted  rectangular  heat  sources  in  three  columns  of  15  elements  per  column. 
The  test  surface  is  immersed  in  a  one  cubic  meter,  plate  glass  tank  filled  with  deionized 
water.  A  parallel  shrouding  plate,  placed  in  front  of  the  test  surface  to  simulate  a  typical 
circuit  board  assembly,  forms  a  channel  with  the  test  surface.  Each  heater  has  a 
copper-constantan  thermocouple  attached  to  its  back.  All  the  thermocouples  are  refer¬ 
enced  individually  to  an  ice  bath  as  seen  in  Figure  2.  Three  thermocouples  monitor  the 
tank  temperature  at  the  top,  middle  and  bottom.  A  Hewlett  Packard  Model  3852S  data 
acquisition  system  measures  the  ice  referenced  voltage  of  each  thermocouple  to  calculate 
temperatures  and  voltages  delivered  from  the  source  and  supplied  to  the  heater  to  cal¬ 
culate  power.  The  voltages  are  converted  into  temperatures  and  component  input 
powers  monitored  by  a  Hewlett  Packard  Model  9153  computer  interfaced  to  the  data 
acquisition  system.  A  filtration  system  consisting  of  four  cartridges  is  used  to  ensure  the 
purity  and  the  resistivity  of  at  least  0.7  megaohm.cm  of  the  water  inside  the  tank. 

A.  MODIFICATIONS  IN  APPARATUS 

The  most  significant  new  feature  was  the  design  of  a  pow'er  assembly  to  provide 
pulsating  input  powers  to  the  heating  elements.  Figure  2  shows  a  sketch  of  this  assem¬ 
bly.  All  the  heater  strips  in  the  left  and  right  columns  are  powered  by 
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SECTION  AA 


Figure  1. 


Test  Assembly  (After  Gaiser  [Ref.  10]).  All  Dimensions  are  in  cm. 
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a  0-40  V,  0-50  A  HP  6269B  D.C.  power  supply  unit  through  four  distribution  panels. 
In  the  central  column  all  heater  strips  except  the  one  at  the  bottom  are  powered  by  a 
0-10  V,  0-10  A  HP  6282A  D.C.  power  supply  unit  through  two  distribution  panels. 
There  is  a  2.01  ohm  precision  resistor  in  series  with  each  heater  resistor  for  current 
measurements.  Details  of  the  power  distribution  assembly  are  provided  by  Gaiser  [Ret. 

10]. 

A  HP-85  computer,  HP-59501B  power  supply  programmer  and  0-40  V,  0-1.5  A 
HP-6289A  D.C.  power  supply  units  were  added  to  the  experimental  setup  to  generate  a 
periodic  power  input.  The  transient  power  input  is  provided  only  to  the  heater  at  the 
bottom  of  the  central  column.  Appropriate  programs  were  developed  for  the  HP-8o 
computer  to  provide  transient  power  to  the  bottom  heater. 

The  following  relation  was  used  to  determine  the  power: 

(As  seen  from  Figure  3) 


POWER  = 


Rp 


X 


where  the  symbols  are  defined  in  the  nomenclature.  See  Appendices  B-F  for  program 
Usiings  for  power  input  and  data  reduction,  developed  for  both  steady  and  transient 

cases. 
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Figure  2.  Experimental  Apparatus 
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Figure  3.  Po»er  Distribution  to  The  Heater  at  The  Bottom  of  Tlie  Centeral  Column 
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B.  EXPERIMENTAL  PROCEDURE 

The  preperation  for  each  set  of  experiments  started  by  stirring  the  water  in  the 
immersion  tank  with  a  propeller  type  stirrer  for  5  to  >0  minutes  to  achieve  uniform 
temperature  distribution.  Crushed  ice  was  placed  into  the  thermocouple  reference  bath. 
A  mercury  thermometer  was  also  inserted  into  the  ice  bath  to  maice  sure  that  the  bath 
was  at  0  ‘’C+0.05  °C  .  The  data  acquisition  system  was  turned  on  and  warmed  up. 
Typically  in  about  one  hour  th:  data  acquision  system  warmed-up  and  also  the 
immersion  tank  was  quiescent.  The  temperature  acquisition  program  (see  Appendix  D) 
was  loaded  and  the  temperatures  of  the  heater  surfaces  and  immersion  tank  were  meas¬ 
ured,  before  the  power  was  supplied  to  the  heaters.  Typically  all  heaters  and  the  tank 
were  within  0.05  °C  of  each  other  prior  to  power  input. 

The  experimental  study  consisted  of  three  major  parts.  The  heater  surface  temper¬ 
ature  measurements  for  three  different  cases  were  made  as  presented  in  Table  1.  In  the 
first  case  thermocouple  measurements  w'ere  taken  every  fifteen  minutes  until  successive 
measurements  remained  unchanged.  Once  steady  state  was  achieved,  temperature 
measurements  and  heater  voltage  readings  were  recorded  and  saved  for  further  process¬ 
ing.  The  computer  program  to  generate  a  prescribed  periodic  power  input,  contumed  in 
Appendix  B  was  then  loaded  and  run.  Once  a  steady  periodic  response  was  evident, 
temperature  oscillations  with  time  were  sampled  over  a  period  of  forty  seconds  with  a 
sampling  frequency  of  180/40  Hz.  The  inpui  pov.cr  pattern  was  also  recorded  over  a 
period  of  seventy  seconds  with  a  sampling  frequency  of  350/70  Hz. 

As  the  amplitudes  of  the  powder  variation  increased,  the  accuracy  in  digitization  of 
the  input  pow'er  by  the  data  acquisition  system  decreased.  While  the  power  input  to  the 
heater  was  as  programmed,  the  measurement  of  it  by  the  HP-85  was  sometimes  slow'er 
than  required.  Another  problem  encountered  was  the  monitoring  of  the  voltage  drops 
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across  both  the  power  supply  and  the  resistor  at  the  same  time  in  order  for  the  calcu¬ 
lation  of  the  power  supplied  to  the  heater. 

In  the  second  catagory  of  experiments  in  Table  1,  the  power  to  all  heaters  was  in¬ 
creased  incrementally.  The  output  from  the  power  calculation  program  was  checked, 
until  the  desired  power  level  was  reached.  The  power  and  temperature  measurements 
were  recorded  when  steady  state  was  attained.  The  periodic  input  power  program  was 
then  loaded  and  run.  Besides  monitoring  the  temperature  oscillations  with  time  on  the 
bottom  heater,  downstream  elements  with  uniform  power  were  also  sampled  to  detect 
any  temperature  oscillations  with  time,  over  a  period  of  five  minutes. 

In  the  third  class  of  experiments  in  Table  1,  the  same  procedure  as  for  the  second 
set  was  repeated.  The  temperatures  of  each  heater  element  both  for  steady  and  transient 
cases  were  recorded.  Additional  runs  for  different  average  power  levels  and  with  differ¬ 
ent  periods  of  applied  power  were  also  performed.  For  the  first  group  of  these  runs,  the 
effects  of  changing  the  pulsation  frequency  and  for  the  second  group,  the  effects  of 
changing  the  amplitude  of  input  power  were  examined.  For  the  third  group  of  runs,  the 
effects  of  providing  all  the  heaters  with  power  different  than  supplid  to  the  bottom  one 
were  studied. 
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Table  1.  EXPERIMENTAL  PROCEDURE  FOR  DIFFERENT  CASES 


CASE 

PROCEDURE 

POWER 

(Watts) 

SHROUD 

SPA.(ntm) 

CASE  1;  Single 

heater  powered  i 

j 

•  Constant  power  is  ap¬ 
plied  to  the  heater  at 
the  bottom  of  centeral 
column,  until  steady 
state  is  achieved. 

•  The  power  input  is  then 
modulated  in  time. 

.049,  .082,  .2, 

.235,  .4,  .6,  .8, 

1. 

Xo  shroud, 

3 

CASE  2:  Centeral 

•  Constant  power  is  ap¬ 
plied  to  all  heaters  until 

.2,  .4,  .6,  .8, 

Same  as 

column  powered 

steady  state  is  achieved. 

•  The  power  to  the  heater 
at  the  bottom  of  cen¬ 
teral  column  is  then 
timewise  modulated 
while  keeping  the  power 
to  all  other  heaters 
steady. 

1.0 

CASE  1 

CASE  3:  All  heat¬ 
ers  powered 

Same  as  CASE  2 

Same  as 

CASE  2 

i 

No  shroud 

(steady  and 

transient), 

3,6,9  (aU 

steady) 
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Table  2.  EXPERIMENTAL  PROCEDURE  FOR  SPECIAL  RUNS 


CASE 

PROCEDURE 

POWER  (Watts) 

Change  in  period  of  in¬ 
put  power 

CASES  1  and  2  in  Table  1 

are  repeated  by  changing  the 

period  from  20  sec.  to  40  sec. 

and  then  10  sec. 

.2,  .6 

Change  in  amplitude 

of  applied  power 

CASE  1  and  2  in  Table  1  are 

repeated  by  changing  the 

amplitude  of  applied  power 

to  selected  values. 

.4,  .8 

Different  power  to 

bottom  heater  than  the 

others  in  central  col¬ 
umn 

Same  as  CASE  2  in  Table  1 

.024  (bot)-  .275  (others), 

.049  (bot)-  .45  (others), 

.082  (bot)-  1.27  (others), 

.235  (bot)-  2.06  (others) 
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III.  EXPERIMENTAL  RESULTS 

A.  SINGLE  HEATER  ELEMENT  DATA 

The  presentation  and  discussion  of  the  single  heater  element  data  during  this  study 
is  divided  into  two  major  parts.  The  first  part  refers  to  the  case  where  the  heater  element 
is  powered  with  uniform  power,  until  steady  state  condition  is  reached.  Within  this  part, 
heater  tempe.ature  data  from  the  thermocouple  at  the  center  of  the  heater  surface  are 
recorded  for  40  sec.  with  a  sampling  frequency  of  180/40  Hz.  In  the  second  set  of  ex¬ 
periments  after  an  initial  steady  state  condition  is  reached  within  the  heater  element,  its 
input  power  is  switched  to  a  sawtooth  form.  In  this  form,  input  power  is  increased  lin¬ 
early  then  decreased  in  the  same  fashion  between  maximum  and  mimumum  amplitudes. 
Applied  transient  power  has  a  period  of  20  sec.  The  surface  temperature  of  the  heater 
element  is  recorded  using  the  the  same  sampling  frequency. 

These  two  sets  of  measurements  are  repeated  for  five  average  power  levels,  ranging 
from  0.2  W  to  I.O  W.  Measurements  are  made  without  shroud  and  also  with  a  shroud 
at  3  mm  spacing  from  the  test  surface,  in  order  to  determine  the  effects  of  power  level 
and  channel  spacing  on  the  heat  transfer  performance.  In  order  to  compare  measured 
surface  temperatures  with  the  numerical  computations,  the  heater  center  temperatures 
are  presented  as  differences  between  the  heater  surface  and  the  ambient. 

1.  Measurements  in  The  Absence  of  Shroud 

In  Figures  4-8  the  temperature  variations  with  time  are  displayed  for  various 
average  power  inputs  with  a  fixed  period  of  20  sec.  In  Figures  9-12  the  effects  of 
changing  the  period  of  applied  power,  for  0.2  W  and  0.6  W  average  power  levels  are 
examined.  In  Figures  13-16  the  effects  of  changing  the  amplitude  of  the  applied  power 
to  0.4  W  and  0.8  W  are  presented.  From  Figures  4-12,  the  trend  of  increasing  average 
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surface  temperature  with  increasing  average  power  is  as  expected.  The  period  of  surface 
temperature  variation,  noted  almost  as  23  sec.  from  these  figures,  stays  approximately 
constant  for  all  power  levels.  It  is  also  seen  that  the  oscillations  in  surface  temperature 
increase  with  an  increase  in  both  amplitude  and  the  period  of  the  applied  power  vari¬ 
ation.  For  example  in  Figures  13  and  14  for  an  average  input  power  of  0.4  \V  when  the 
power  variation  is  within  0.55  W  and  0.29  W,  the  maximum  difference  in  AT  is  O.S  °C. 
This  difference  decreases  to  0.25  °C,  if  the  power  variation  is  kept  within  0.43  W  and 
0.36  W.  For  0.6  W  average  power  level  in  Figures  15-16  when  the  period  of  applied 
power  is  increased  to  40  sec.,  the  difference  in  AT  increases  to  1.0  °C.  The  difference  in 
AT  decreases  to  0.3  °C  when  period  is  reduced  to  10  sec.  For  higher  power  levels,  the 
effect  of  amplitude  and  period  of  applied  power  is  more  evident. 

2.  Measurements  With  Shroud 

Figures  17  and  18  show  AT  versus  time  at  mean  power  levels  of  0.2  W  and  1.0 
W  for  shroud  spacing  of  3  mm.  When  compared  to  the  unshrouded  measurements  at 
the  same  power  levels  [Figures  4  and  8),  it  is  seen  that  the  AT  values  for  a  single  element 
remain  almost  unchanged. 
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TIME  VS  POWER  (tJO  SH) 


TIME  VS  DELTR  T  (NO  SH ) 


Figure  4.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  at  0.2  W.  Pulsation  Amplitude  and  Period  are  0.06  W 
and  20  sec 
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TIME  VS  POWER  (NO  SH) 


Figure  6 


TIME  VS  DELTR  T  (NO  SH) 


Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  at  0.6  W.  Pulsation  Amplitude  and  Period  are  0.1  \V 
and  20  sec. 


BHrrna  h  h  3)  s:  xd  tn  2:  o -□ 


TIME  VS  POWER  (NO  SH) 


Figure  8 


Power  Pulsations  at  1.0  W.  Pulsation  Amplitude  and  Period  are  0.13  W 
and  20  sec 
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Figure  9, 


TIME  VS  POWER  (NO  SH ) 


time  vs  DELTH  T  (NO  SH) 


Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 

Power  Pulsations  with  a  Period  of  10  sec.  at  0.2  W.  Pulsation  Amplitude 
is  0.06  W. 
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Figure  11 


TIME  VS  POWER  (NO  SH ) 


TIME  VS  DELTFI  T  (NO  SH) 


Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Po^ver  Pulsations  with  a  Period  of  10  sec.  at  0.6  W.  Pulsation  Amplitude 
is  0.1  W. 
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Figure  12.  Oscilla 
Power  P 


is  0.1  W 


TIME  VS  POWER  (NO  SH) 


TIME  VS  DELTA  T  (NO  SH ) 


Figure  13.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  with  Mean  Power  of  0.4  W,  Amplitude  of  0.03  W  and 
Period  of  20  sec. 
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Q.  o  2  u  q:  2  tr  I- h-  auiJHo: 


TIME  VS  POWER  (NO  SH ) 


TIME  VS  ]]ELTfJ  I  (NO  SID 


Figure  14.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  with  Mean  Power  of  0.4  W,  Amplitude  of  0.13  W  and 
Period  of  20  sec. 


TIME  VS  POWER  (NO  SH) 


TIME  VS  DELTR  T  (NO  SH) 


Figure  15.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  with  Mean  Power  of  0.8  W,  Amplitude  of  0.06  W  and 
Period  of  20  sec. 
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TIME  VS  POWER  (NO  SH) 


TIME  VS  DELTA  T  (NO  SH) 


Figure  16.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  'nput 
Power  Pulsations  with  Mean  Power  of  0.8  W,  Amplitude  of  0.19  W  and 
Period  of  20  sec. 
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TIME  VS  DELTR  T  (3  MM  SH.  SP . ) 


Figure  17.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  with  Mean  Power  of  0.2  W,  Amplitude  of  0.06  W  and 


Period  of  20  sec. 

TIME  VS  DELTA  T  (3  MM  SH.  SP . ) 


Figure  18.  Oscillations  in  Surface  Temperature  of  Bottom  Heater  Due  To  Input 
Power  Pulsations  with  Mean  Power  of  1.0  W,  Amplitude  of  0.13  W  and 
Period  of  20  sec. 
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B.  SINGLE  COLUMN  DATA 

In  another  set  of  experiments  the  efiects  of  transient  power,  applied  to  the  heater 
at  the  bottom  of  the  central  column  were  examined.  The  remaining  fourteen  heaters  in 
the  column  were  steadily  powered.  The  main  point  of  interest  was  to  look  for  changes 
in  the  heat  transfer  characteristics  compared  with  those  existing  with  steady  power  to 
all  heaters.  The  same  experimental  procedure  that  is  followed  in  the  previous  section  is 
used  for  powering  the  bottom  heater.  The  steady  and  transient  temperature  data  were 
taken  for  five  power  levels,  ranging  from  0.2  W  to  1.0  W  both  with  and  without  a 
shroud.  The  surface  temperatures  of  heaters  1st,  10th  and  13th  were  monitored  once 
steady  power  of  bottom  heater  w'as  switched  to  the  periodic  form.  Temperature  meas¬ 
urements  were  made  with  a  sampling  frequency  of  300/335  Hz.  over  a  period  of  335 
seconds. 

1.  Measurements  In  The  Absence  of  Shroud 

Figures  19-23  show  AT  versus  x/L  at  steady  input  power  levels  of  0.2  W,  0.4 
W,  0.6  W,  0.8  W  and  1.0  W  respectively.  From  the  plots  of  data,  the  trend  of  increasing 
surface  temperature  along  the  test  surface  is  readily  apparent.  It  is  also  seen  that  at  a 
given  location,  surface  temperature  increases  with  rising  power  level. 

The  steady  temperature  levels  were  used  to  determine  the  non-dimensional  heat 
transfer  characteristics.  Rayleigh  and  Nusselt  numbers  were  calculated  as; 

Ra  =  igPifP)  /  (akfV),  Nu  =  h  \  /  k^.  Figures  24  and  25  are  plots  of  Nu  versus  Ra 
for  0.2  W  and  0.8  W  respectively.  The  local  values  of  h  indicate  a  decreasing  trend 
moving  from  the  bottom  heater  to  the  top.  It  is  also  obvious  that  h  values  increase  with 
rising  powder  level  at  a  given  location  due  to  the  stronger  flow. 

The  transient  temperature  data  for  the  bottom  heater  over  a  period  of  40  sec. 
were  the  same  as  for  the  single  heater  element  examined  in  the  previous  section.  Figures 
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26  and  27  are  the  temperatures  of  1st,  lOth  and  I3th  heaters,  over  a  period  of  335  sec. 
for  0.4  W  and  0.8  VV  input  power  levels  respectively.  As  seen  from  the  figures,  oscillation 
in  surface  temperatures  of  downstream  heaters,  when  power  to  the  bottom  heater  is  os¬ 
cillated  with  time  are  barely  noticeable. 

2.  Measurements  With  Shroud 

In  the  previous  section,  it  was  found  that  the  effects  of  the  shroud  on  the  heat 
transfer  from  the  single  heater  element  were  not  apparent.  In  this  section,  the  effect  of 
shroud  is  investigated  with  the  entire  central  column  powered.  Surface  temperature 
measurements  of  heaters  were  recorded  at  five  different  power  levels  for  3  mm  shroud 
spacing.  The  results  are  compared  with  the  unshrouded  measurements. 

The  temperature  rise  above  ambient  for  each  heater  at  average  power  levels  of 
0.2  W  and  0.8  W  for  3  mm  shroud  spacing  are  seen  in  Figures  28  and  29.  From  these 
plots  temperature  values  higher  than  the  unshrouded  case  measurements  are  seen  be¬ 
yond  the  sixth  heater  element  from  the  bottom.  When  the  power  levels  increase,  difi'er- 
ences  in  AT  values  also  increase  compared  to  the  unshrouded  case.  These  measurements 
are  in  agreement  with  Knight  (Ref  13]  and  Gaiser  [Ref  10).  While  the  thermal  layer 
near  the  bottom  component  is  thin  enough  not  to  be  effected  by  the  presence  of  the 
shroud,  the  same  is  not  true  further  downstream. 
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26  and  27  are  the  temperatures  of  1st,  10th  and  13th  heaters,  over  a  period  of  335  sec. 
for  0.4  W  and  0.8  W  input  power  levels  respectively.  As  seen  from  the  figures,  oscillation 
in  surface  temperatures  of  downstream  heaters,  when  power  to  the  bottom  heater  is  os¬ 
cillated  with  time  arc  barely  noticeable. 

2.  Measurements  With  Shroud 

In  the  previous  section,  it  was  found  that  the  eflecis  of  the  shroud  on  the  heat 
transfer  from  the  single  heater  element  were  not  apparent.  In  this  section,  the  en'ect  of 
shroud  is  investigated  with  the  entire  central  column  powered.  Surface  temperature 
measurements  of  heaters  were  recorded  at  five  different  power  levels  for  3  mm  shroud 
spacing.  The  results  are  compared  with  the  unshrouded  measurements. 

The  temperature  rise  above  ambient  for  each  heater  at  average  power  levels  of 
0.2  W  and  0.8  W  for  3  mm  shroud  spacing  are  seen  in  Figures  28  and  29.  From  these 
plots  temperature  values  higher  than  the  unshrouded  case  measurements  are  seen  be¬ 
yond  the  sixth  heater  element  from  the  bottom.  When  the  power  levels  increase,  difler- 
ences  in  AT  values  also  increase  compared  to  the  unshrouded  case.  These  measurements 
are  in  agreement  with  Knight  (Ref.  13)  and  Gaiser  (Ref  10).  While  the  thermal  layer 
near  the  bottom  component  is  thin  enough  not  to  be  effected  by  the  presence  of  the 
shroud,  the  same  is  not  true  further  downstream. 
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X/L  VS  DELTA  T  (C) 


XL"L  VS  DELTA  T  (C) 

11. 1 


X^L 

Figure  22.  x/L  versus  AT  at  0.8  W  Input  Power 
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RR  VS  NU 


Figure  24.  Nu  versus  Ra  at  0.2  W 


RR  VS  tJ’J 


RR 

Figure  25.  Nu  versus  Ra  at  0.8  W 
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TIME  VS  DELTn  T  (C)  (NO  SH) 


Figure  26. 


Figure  27. 


Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.4  W 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.07  W  and  20  sec.  respectively. 

TIME  VS  DELTR  T  (C)  (NO  SH) 
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Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.8  W 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.07  W  and  20  sec.  respectively. 
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C.  DATA  WITH  ALL  THREE  COLUMNS  POWERED 


A  third  set  of  experiments  was  undertaken  to  examine  the  effects  of  placing  flanking 
columns  of  heaters  on  each  side  of  the  central  column.  The  three  column  results  are 
discussed  next.  In  addition  to  the  measurements  for  the  central  column,  the  surface 
temperatures  of  top  and  bottom  heaters  of  both  side  columns  were  also  monitored  for 
335  sec.  with  a  sampling  frequency  of  300  335  Hz.  These  transient  measurements  were 
repeated  for  average  power  levels  ranging  from  0.2  W  to  1 .0  W  with  3mm  shroud  spacing 
and  also  without  shroud.  Additional  runs  were  performed  with  steady  power  inputs, 
within  the  same  power  range  and  shroud  spacings  of  6  and  9  mm. 

1.  Measurement.s  In  The  Absence  of  Shroud 

The  surface  temperatures  in  each  of  the  three  columns  are  plotted  to  indicate 
changes  in  convective  heat  transfer  characteristics  as  compared  with  the  single  column 
case.  Figures  .’0  and  31  are  AT  versus  x/L  plots  at  steady  power  levels  of  0.2  W  and  0.8 
W  respectively,  from  these  plots  it  is  seen  that  at  any  x/L  the  central  column  has 
slightly  higher  AT  values  when  compared  to  the  right  and  left  columns.  Referring  back 
to  Figures  19  and  23  for  a  single  column,  it  is  apparent  that  the  AT  values  for  central 
column  components  with  all  three  columns  powered  are  similar  to  the  values  obtained 
for  single  column  powered  case.  Comparing  the  heat  transfer  coefficients  for  these  ex¬ 
periments  with  the  single  column  values  it  was  seen  that  the  three  column  data  produced 
almost  the  same  values  of  h. 

As  was  done  with  the  single  column,  the  surface  temperature  of  the  bottom 
heater  of  the  central  column  was  monitored  for  40  sec.  with  a  sampling  frequency  of 
180/40  Hz,  after  modulating  its  input  power  with  a  time  period  of  20  sec.  0’^«'e  again 
the  temperatures  of  the  three  heaters  above  it  and  the  heater  at  the  top  of  each  column 
were  also  monitored  for  335  sec.  with  a  sampling  frequency  of  300/335  Hz.  Figures 
32-34  are  plots  of  AT  versus  time  for  mean  power  0.2  W,  amplitude  0.06  W  and  Figures 
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35-37  for  average  power  0.8  W,  amplitude  0.12  W.  From  the  figures,  no  oscillation  in 
the  surface  temperature  in  any  heater  indicates  the  damping  of  the  effect  of  the  input 
power  variation  in  the  bottom  heater  is  observed.  It  is  also  noted  that  the  surface  tem¬ 
perature  oscillation  of  bottom  heater  due  to  periodic  power,  follows  the  same  pattern 
as  for  a  single  heater. 

2.  Measurements  With  Shroud 

-Measured  heater  surface  temperatures  are  presented  for  the  central  column  in 
Figures  38-40  at  various  shroud  spacings  and  power  levels.  It  is  seen  that  the  shroud 
has  a  significant  effect  on  the  heater  temperatures  only  for  small  spacings.  For  spacings 
greater  than  6  mm,  the  heater  temperatures  are  almost  equal  to  the  unshrouded  test 
case.  Figure  41  compares  the  surface  temperatures  for  3  mm  shroud  spacing  at  various 
power  levels.  Figure  42  compares  the  measured  surface  temperatures  without  shroud 
with  the  values  for  3  mm  shroud  spacing  at  0.4  W  and  0.8  W  power  levels.  From  the 
figure  it  is  seen  that  as  the  shroud  spacing  decreases,  surface  temperatures  of  central 
column  heaters  are  increased.  Figure  43  shows  the  surface  temperatures  of  all  heaters 
for  3  mm  shroud  spacing  at  0.4  W.  From  the  figure  it  is  seen  that  for  this  shroud  spac¬ 
ing,  surface  temperatures  of  the  central  column  heaters  are  higher  than  the  measured 
values  for  the  side  columns  at  the  same  x,'L.  Fluid  entrainment  with  shroud  in  place  is 
impeded  in  the  direction  normal  to  the  test  surface.  Lateral  fluid  movement  brings 
warmer  fluid  from  the  flanking  columns  towards  the  central  column.  This  results  in  an 
increase  in  temperature  of  the  central  column  heaters. 
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TIME  VS  DELTR  T  (C)  (fJO  SH) 
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Figure  32.  Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.2  W 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.06  W  and  20  sec.  respectively. 


TIME-SEC 


Figure  33.  Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.2  W 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.06  W  and  20  sec.  respectively. 


TIME  VS  DELTR  T  (C)  (NO  SH) 
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Figure  34.  Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.2  W 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Perio..  are  0.06  W  and  20  sec.  respectively. 


TIME  VS  DELTA  T  (C)  (NO  SH) 


Figure  35.  Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.8  W 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.12  W  and  20  sec.  respectively. 


TIME  VS  DELTA  T  (C)  (NO  SH ) 


Figure  36.  Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.8  W' 
Wlien  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.12  W  and  20  sec.  respectively. 
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TIME  VS  DELTA  T  (C)  (NO  SH) 


TIME-SEC 

Figure  37.  Monitoring  The  Surface  Temperatures  of  Selected  Heaters  at  0.8  VV 
When  Input  Power  To  The  Bottom  Heater  Oscillated.  Oscillation  Am¬ 
plitude  and  Period  are  0.12  W  and  20  sec.  respectively. 
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Figure  38.  x/L  versus  AT  for  Different  Shroud  Spacings  At  0.2  W 


44 


X/L  VS  DELTA  T 


> 


LEnESD 
o  3  MM 

'c  r.  MSI 

i  9  MM 


POWER  =  0.6  WATTS 


0.0  0  1  0.2  0.3  0  4  0  3 

X/L 


Figure  39.  x/L  versus  AT  for  Different  Shroud  Spacings  At  0.6  VV 


X/L  VS  DELTA  T 


Figure  40.  x/L  versus  AT  for  Different  Shroud  Spacings  At  1.0  W 
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X/L  VS  DELTA  T 


Figure  41.  x/L  versus  AT  for  3  mm  Shroud  Spacing  At  Various  Power  Levels 


X/L  VS  DELTA  T 


Figure  42.  Comparison  of  Surface  Temperatures  Measured  Without  And  With  A 
Shroud  at  3  mm  For  0.4  W  and  0.8  W 
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IV.  NUMERICAL  MODEL  AND  COMPARISON  WITH 

MEASUREMENTS 

Based  on  the  experimental  results  presented  in  the  previous  chapter  it  is  clear  that 
in  the  absence  of  shroud  the  lateral  fluid  entrainment  is  not  ver>’  strong.  With  only  a 
single  component  powered,  the  transport  away  from  the  ends  can  then  be  modelled  as 
two  dimensional.  In  the  following,  a  computational  model  of  natural  convection  trans¬ 
port  from  a  single  flush  heater  is  presented. 

A.  MATHEMATICAL  FORMULATION 

The  schematic  sketch  of  the  configuration  numerically  investigated  is  showm  in  Fig¬ 
ure  44.  The  flow  is  assumed  to  be  two  dimensional,  resulting  from  a  heat  source,  flush 
mounted  on  a  vertical  substrate  and  located  at  the  middle  of  the  surface.  The  shaded 
region  represents  the  substrate,  and  the  darker  region  in  it  represents  the  heater.  The 
heater,  the  substrate  and  the  fluid  have  constant  thermophysical  properties.  The  vertical 
coordinate  x  is  measured  upw'ard  from  the  base  of  substrate  and  u  is  the  velocity  com¬ 
ponent  in  this  direction.  The  y  coordinate  is  measured  horizontally  from  the  surface  and 
V  is  the  velocity  component  in  the  same  direction.  Isothermal  impermeable  enclosure 
walls  with  no-slip  i.e  u  =  v  =  T  =  0  are  prescribed  as  boundary  conditions.  The  governing 
equations  assuming  no  viscous  dissipation  and  pressure  energy  terms  and  the  Boussinesq 
approximations  to  be  true,  are: 

FLUID  REGION : 


48 


5(p/)  G{pfU)  d{pfV) 

+  - :r -  +  - :r -  =  0 


0  t 


Sy 


d{pfu)  c{pfu)  c{pfUv) 

+  - r -  + 


c  [ 


C  X 


oy 


,  u  .  u  ^ 
=  ) 


O  X 


oy 


+  pfg  P{T  -  Tamb  )  - 


d{pfV)  d{p^uv)  d{pfV  ) 

+  - -  +  - ; -  =  p{ 


o  l 


0  X 


cy 


V  ^  ^  P 

d  x^  d  ^  y 


d{pfCpT)  d{pfCpuT)  d{pfCpvT)  .  ,  T  T 


dt 


+ 


d  X 


^y 


+ 


C  X 


SUBSTRATE: 


^iPs<^pT)  ,  ,  T  T 

“TT-  "*■(77  ^  77 


HEATER  REGION  : 


diPfjCp  T) 
dt 


( 


ilZ 

dx^ 


+ 


d^  T 


+  <1 


/// 


49 


B.  COMPUTATIONAL  PROCEDURE 


The  information  contained  in  the  differential  equations  is  discretized  which  makes 
it  possible  to  replace  the  governing  partial  differential  equations  with  algebraic 
equations,  which  can  be  solved  with  relative  ease.  A  finite  difference  scheme  is  used  for 
discretizing  the  equations.  The  calculation  domain  is  divided  into  a  number  of  control 
volumes  such  that  there  is  one  control  volume  surrounding  each  grid  point.  The  tem¬ 
perature  and  pressure  components  are  located  at  the  main  grid  points.  The  location  of 
the  velocity  components  is  staggered  and  defined  on  the  control  volume  faces.  The 
power  law  profile  that  provides  a  ver>'  good  fit  to  the  exact  expression,  is  used  for  the 
spatial  variation  of  the  dependent  variables  between  grid  points.  The  harmonic  mean 
formulation  that  provides  physically  meaningful  results  even  with  large  changes  in  ma¬ 
terial  properties  is  followed  for  the  calculation  of  the  interface  diffusivities.  The  details 
of  the  discretization  process  can  be  found  in  Patankar  [Ref.  14). 

The  discretized,  nominally  linear  equations  are  solved  by  line-by-line  TDMA  (  Tri- 
Diagonal  Matrix  Algorithm  )  method  in  which  the  coefficients  are  iteratively  recalcu¬ 
lated  and  the  process  repeated  until  convergence.  ‘  Convergence  is  based  on  a  balance 
of  the  rate  of  energy  generated  in  the  heater  and  the  rate  of  energy  leaving  the  enclosure 
walls.  The  SIMPLER  (SIMPLE- Revised)  algorithm  that  employs  two  different 
equations:  one  for  the  task  of  correcting  the  velocities,  the  other  for  obtaining  a  new 
pressure  field,  is  followed  during  solution  as  outlined  in  [Ref.  14).  Even  though  separate 
equations  are  written  for  solid  and  fluid  regions,  numerically  the  solid  region  in  the  mo¬ 
mentum  equation  is  simulated  by  specifying  a  very  large  viscosity.  The  heat  source  term 
and  property  values  in  the  energy  equation  are  specified  appropriately  in  the  various 
regions.  As  a  result,  only  one  momentum  and  one  energy  equation  are  solved  through¬ 
out  the  computational  domain. 
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Numerical  computations  were  performed  on  48x56  control  volumes.  Different  time 
steps  were  tried  ranging  from  0.1  sec.  to  5.0  sec.  to  determine  their  effect  on  the  sol¬ 
utions.  Calculations  for  the  modulated  power  input  were  started  with  u,  v  and  T  values, 
determined  previously  for  a  steady  case  with  the  same  average  power.  The  initial  profiles 
were  obtained  by  running  the  transient  code  with  a  very  large  time  step. 

The  computations  were  performed  initially  on  a  IBM-370;3033  computer  and  later 
on  an  Amdahl  Computer.  Appendix  1  shows  the  problem  dependent  portion  of  the  fi¬ 
nite  difference  program  used  to  obtain  the  the  numerical  results. 
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C.  COMPARISON  OF  COMPUTATIONS  WITH  EXPERIMENTAL 


MEASUREMENTS 

Numerical  investigation  of  transient  heat  transfer  consisted  of  two  parts.  For  the 
first,  initial  transient  response  due  to  sudden  step  power  input  was  studied.  Figures  45 
and  46  show  representative  transient  How  and  temperature  contour  developments  for 
8.69  W  m  step  input.  We  note  that  for  the  heater  size  used  in  the  experiments  this 
power  input  corresponds  to  0.2  W.  The  figures  from  left  to  right  show  the  contours  at 
10,  44,  so  and  120  seconds  following  start  of  heating.  At  10  sec.  the  isotherms  are  al¬ 
most  parellel  to  the  heater  indicating  a  diffusion  dominated  response.  A  single  clockwise 
cell  exists  next  to  the  heater.  In  the  second  and  third  figures,  the  secondarx'  cell  appears 
and  grows.  Flow  becomes  relatively  faster  and  spreads  further  away.  Also  at  120  sec. 
the  flow  begins  to  resemble  the  fully  developed  condition  seen  in  Figure  54  for  the  same 
input  power. 

For  the  second  part  of  numerical  study,  the  transient  response  following  a  periodic 
input  power  into  the  heater  was  studied.  The  amplitude  and  time  period  of  imposed 
power  oscillations  were  varied.  Time  steps  from  0.1  sec.  to  5.0  sec.  were  investigated. 
Surface  temperatures  calculated  for  different  time  steps  are  presented  in  Table  3.  Tran¬ 
sient  variations  of  surface  temperatures  for  different  time  steps  are  displayed  in  Figures 
47  -  50.  The  surface  temperature  patterns  are  similar  for  all  time  steps.  The  highest 
difference  in  calculated  temperature  values  was  2.25  %,  when  5.0  sec.  time  step  was  used 
instead  of  0. 1  sec.  Difference  in  calculated  surface  temperatures  was  smaller  than  0.45 
°/'o,  if  calculations  were  performed  with  0.5  sec.  time  step  rather  than  0.1  sec.  In  order 
to  optimize  the  computation  times  without  compromising  significantly  on  the  accuracy, 
all  numerical  calculations  were  carried  out  with  a  time  step  of  0.5  sec. 

Numerical  computations  were  also  performed  for  two  different  sets  of  heating  ele¬ 
ment  properties  since  the  actual  element  properties  in  the  experiments  were  not  easily 
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quantifiable.  In  the  experiment,  the  composite  heater  element  consisted  of  an  inconel 
ribbon  attached  to  a  kapton  substrate.  In  order  to  cover  the  range  of  property  vari¬ 
ations  between  the  two  materials,  the  computations  were  repeated,  with  the  properties 
of  inconel  and  of  plastic  separately  for  the  heater  [Ref  15].  Figures  51  and  52  show  the 
numerically  calculated  surface  temperatures  with  time  for  average  power  levels  of  0.2 
W  and  1.0  W  for  plastic  heater  and  Figure  53  for  inconel  heater  at  1.0  W  average  power 
level.  While  the  mean  temperatures  in  Figure  53  are  lower  than  in  Figure  52,  the  am¬ 
plitudes  of  temperature  oscillations  remain  approximately  the  same. 


Table  3.  EFFECT  OF  SELECTING  DIFFERENT  TIME  STEPS  ON  COM¬ 


PUTED  QUANTITIES 


TIMEJsec.) 

Surface  Temperatures  (°Q  Calculated  by  Taking  Time  Steps  (sec.)  of  ; 

0.1 

0.5 

1.0 

5.0 

0.1 

2.1271 

- 

- 

- 

0.5 

2.1186 

2.1179 

- 

- 

1.0 

2.1094 

2.1081 

2.1072 

- 

2.0 

2.0955 

2.0930 

2.0910 

- 

5.0 

2.0866 

2.0806 

2.0756 

2.0394 

10.0 

2.1615 

2.1512 

2.1412 

2.0660 

15.0 

2.2149 

2.2104 

2.2074 

2.1848 

20.0 

2.1474 

2.1484 

2.1522 

2.1802 

Table  4.  THERMOPHYSICAL  PROPERTIES  OF  MATERIALS  IN  HEATER 


ELEMENT 


Thermophysical 

Property 

MATERIAL 

Inconel 

Plastic 

c,(J/kg-K) 

439.0 

1500.0 

k  (W;m-K) 

11.7 

0.2 

P  (kg/m’) 

8510.0 

1180.0 
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Figure  45.  Development  of  Flow  During  Initial  Transient  Due  To  0.2  W  Step  Input 
Power.  The  Figures  are  at  (i)10,  (ii)44,  (iii)80  and  (iv)120  seconds. 
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SURfACK  TEMP  vs  TIMB 
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Figure  47.  Surface  Temperature  Calculated  with  Time  Step  of  0. 1  sec. 
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SURFACE  TEMP  VS  TIME 
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Figure  50.  Surface  Temperature  Calculated  with  Time  Step  of  5.0  sec. 


0.0  9  0  10.0  15.0  20.0 


riME-(sf:c) 

Figure  51.  Surface  Temperature  Calculated  with  the  Properties  of  Plastic  at  0.2  W 


Average  Input  Power 
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Figure  52.  Surface  Temperature  Calculated  with  the  Properties  of  Plastic  at  1.0  W 
Average  Input  Power 
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Figure  53.  Surface  Temperature  Calculated  with  Ihe  Properties  of  Inconel  at  1.0 
W  Average  Input  Power 


COMPARISON 

The  bottom  heater  in  the  central  column  (Section  AA  in  Fig.4)  was  chosen  for  a 
comparison  with  numerical  results.  For  0.2  W  average  power  input  with  0.06  W  ampli¬ 
tude  and  20  sec.  pulsations,  the  steady  periodic  temuerature  response,  calculated  at  the 
back  of  the  model  heater,  was  27  %  lower  than  the  measured  for  an  inconel  heater  and 
38  %  higher  than  the  measured  value  for  a  plastic  heater.  The  oscillation  amplitudes  in 
calculated  surface  temperatures  agreed  within  approximately  33  %  with  the  exper¬ 
imentally  measured  values.  For  the  1.0  W  average  power  case,  the  calculated  surface 
temperature  of  an  inconel  heater  was  13  %  lower  than  the  experimentally  measured 
value.  The  calculated  value  for  a  plastic  heater  was  51  %  higher  than  the  measured  one. 
The  amplitudes  of  oscillations  in  calculated  surface  temperatures  agreed  within  33  % 
with  the  experimental  values. 

The  percent  difierences,  presented  in  the  following  between  numerically  calculated 
and  experimentally  measured  surface  temperatures,  arc  for  the  case  when  the  heater 
material  is  assumed  to  be  inconel.  Figures  54  and  55  show  the  corresponding  temper¬ 
ature  contours  and  streamlines  for  the  average  power  input  levels  of  0.2  W  and  1.  W. 
The  effect  of  change  in  the  period  of  the  input  power  pattern  was  also  examined  nu¬ 
merically.  Figures  56  and  57  show  the  surface  temperature  variation  with  time  for  0.2 
W  average  power  with  different  periods.  The  trend  for  increasing  period  of  the  measured 
temperature  response  due  to  increase  in  period  of  applied  power  is  evident.  When 
compared  to  the  experimentally  measured  values,  the  calculated  temperature  oscillation 
amplitudes  were  almost  33  %  lower  for  all  pow'er  levels.  Figures  58-60  show  the  effects 
of  the  period  of  applied  power  for  a  higher  average  power  level  of  0.6  W. 

The  effect  of  change  in  the  amplitude  of  the  applied  power  variation  was  also  in¬ 
vestigated.  Computations  revealed  increasing  differences  between  peak  values  of  calcu- 
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lated  heater  surface  temperatures  due  to  increase  in  amplitude  of  applied  power.  The 
calculated  surface  temperatures  versus  time  for  0.4  W  and  0.8  W  average  power  levels 
with  changing  amplitudes  are  presented  in  Figures  61  and  62.  The  corresponding  tem¬ 
perature  contours  and  streamlines  are  plotted  in  Figures  63  and  64.  The  calculated 
steady  mean  temperature  was  almost  21  %  lower  for  0.4  W  and  13  %  lower  for  0.8  W 
compared  to  the  experimentally  measured  values.  The  calculated  oscillation  amplitudes 
in  surface  temperatures  agreed  only  within  67  %  with  the  experimental  values  for  both 
power  levels. 

This  difference  can  be  ascribed  to  the  larger  thickness  of  the  heater  model  than  the 
actual  value.  The  actual  thickness  was  about  .195  mm.  However  a  very'  fme  grid  will 
be  required  if  this  size  is  to  be  taken  in  the  model.  As  the  thickness  of  the  heater  in  the 
computational  model  was  reduced,  the  oscillation  amplitudes  increase.  This  trend  is 
presented  in  Figure  66  with  the  depth  of  the  heater  reduced  to  half  of  the  previous  value. 
If  the  actual  properties  of  heater  element  which  are  combination  of  inconel  and  kapton 
properties  were  available,  it  is  obvious  that  the  agreement  with  the  experiment  will  im¬ 
prove,  The  experimentally  measured  temperature  values  are  bracketed  by  the  compu¬ 
tations,  as  expected. 

Additional  computations  were  performed  with  k  values  between  inconel  and  plastic 
to  achieve  better  agreement  between  the  model  and  experiments.  Figure  65  shows  the 
surface  temperatures  calculated  with  k„  =  1.45  It  is  seen  that  calculated  val¬ 

ues  are  reasonably  close  to  the  experimental  values.  A  particular  run  was  performed 
with  k=  1.0  W/m-K  and  heater  thickness  of  .195  mm.  Computed  results  are  presented 
in  Figure  67.  For  these  choices  of  thermophysical  properties  and  component  dimensions 
close  agreement  is  seen  with  the  experimental  data  in  Figure  4. 
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Figure  56.  Surface  Temperature  Due  to  Power  Pulsation  with  a  Period  of  40  sec. 
for  0.2  W  Average  Power  and  Amplitude  of  0.06  W. 
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SURFACE  TEMP  VS  TIME 
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Figure  58.  Surface  Temperature  Due  to  Power  Pulsation  with  a  Period  of  40  sec 
for  0.6  W  Average  Pow  er  and  Amplitude  of  0. 1  W 
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Figure  59.  Surface  Temperature  Due  to  Power  Pulsation  with  a  Period  of  20  sec 


for  0.6  W  Average  Power  and  Amplitude  of  0.1  \V 
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Figure  60.  Surface  Temperature  Due  to  Power  Pulsation  with  a  Period  of  10  sec. 
for  0.6  W  Average  Power  and  Amplitude  of  0.1  W 
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Figure  61.  Surface  Temperature  Due  to  Power  Pulsation  with  Different  Amplitudes 


for  0.4  W  Average  Power 
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SURFACE  TEMP  VS  TIME 
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Figure  62.  Surface  Temperature  Due  to  Po..er  Pulsatiou  with  Different  Amplitudes 


for  0.8  W  Average  Power 


70 


Figure  63.  Temperature  Contours  and  Streamlines  for  0.4  W  Average  Power,  0.07 


W  Amplitude  and  20  sec.  Period  of  Oscillation 
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Figure  65. 


Figure  66. 
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Effects  of  Thermophysical  Properties  on  Calculated  Temperatures 


SURFACE  TEMP  VS  TIME 


TJME-(SrC) 


Effects  of  Heater  Dimensions  on  Calculated  Temperatures 
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Figure  67.  Surface  Temperature  vs  Time 
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V.  CONCLUSIONS 


Following  conclusions  are  drawn  based  on  this  experimental  and  numerical  study 
about  the  heat  transfer  effects  of  input  power  pulsations. 

The  experiments: 

The  oscillation  amplitudes  in  surface  temperature  increase  with  an  increase  in  both 
the  amplitude  and  the  time  period  of  the  applied  power.  The  period  of  temperature  os¬ 
cillation  remains  almost  the  same  for  all  shroud  spacings  at  a  particular  power  level.  It 
was  observed  that  the  effects  of  periodic  power,  applied  to  the  heater  at  the  bottom  of 
the  central  column  are  barely  noticeable  on  the  heat  transfer  characteristics  of  the  other 
steadily  powered  heaters  in  the  column.  There  is  no  noticeable  oscillation  in  surface 
temperatures  of  dowmstream  heaters  w’hen  power  to  the  bottom  heater  oscillates  in  time. 
With  all  three  columns  powered  under  steady  conditions  and  with  a  shroud  in  place,  the 
central  column  heaters  at  any  x/L  have  slightly  higher  AT  values  w’hen  compared  to  the 
right  and  left  columns.  For  spacings  greater  than  6  mm,  the  heater  temperatures  are 
almost  equal  to  the  unshrouded  test  case.  These  conclusions  regarding  the  steady  re¬ 
sponse  are  in  agreement  with  Knight  [Ref.  13)  and  Gaiser  [Ref.  10). 

The  computations: 

Numerical  computations  of  the  transient  response  due  to  both  a  sudden  step  power 
input  and  a  periodic  power  input  were  conducted.  Transient  flow  and  temperature  pro¬ 
files  were  obtained.  Computations  were  performed  for  two  different  heating  element 
properties  since  the  actual  element  in  the  experiments  was  a  composite  and  its  properties 
were  not  easily  quantifiable.  Thermophysical  properties  of  the  heating  element  and  its 
thickness  had  a  pronounced  effect  on  calculated  temperatures.  As  the  thickness  of  the 
heater  element  was  reduced,  the  temperature  oscillation  amplitudes  increased.  Once  the 
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variations  in  the  thermophysical  properties  and  appropriate  thickness  of  the  heating  el¬ 
ement  were  taken  into  account  in  the  numerical  study,  the  predicted  temperatures  were 
in  reasonable  agreement  with  measurements. 
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VI.  RECOMMENDATIONS 

In  continuation  of  this  study  it  is  reconrunended  that  the  following  aspects  be  further 
investigated: 

•  The  eflects  of  pulsating  power  inputs  to  more  than  one  heater  on  the  heat  transfer 
from  the  other  steadily  powered  heaters  should  be  considered. 

•  Investigate  heat  transfer  enhancement  due  to  possible  turbulence  effects  caused  by 
high  amplitude  power  pulsations  at  relatively  high  average  power  levels. 

•  Use  a  multiple  heater  numerical  model  to  determine  heater  interaction  effects. 

•  Analyze  transient  conditions  for  both  one  and  multiple  heater  numerical  models  in 
the  presence  of  a  shroud. 

•  Numerically  study  the  effects  of  thermophysical  properties  of  substrate  on  the  am¬ 
plitudes  of  oscillations  of  surface  temperature. 
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APPENDIX  A.  SUBROUTINE  USER 


SUBROUTINE  USER 

LOGICAL  LSOLVE,LPRINT,LBLK,LSTOP.LSWICH 

COMMON  F(66,B6,10),P(66.66),RHO(66,66),GAM(66.66),CON(66,66), 

1  AIP(66,66).AIM(66,66),AJP(66,66)>AJM(66,66),AP(66,66), 

2  X(66),XU(66),XDIF(66),XCV(66),XCVS(66), 

3  Y(66),YV(66),YDIF(66),YCV(66),YCVS(66), 

A  YCVR(66),YCVRS<66),ARX(66).ARXJ(66),ARXJP(66), 

5  R(66),RMN(66),SX(66),SXMN(66),XCVI(66),XCVIP(66) 

COMMON  DU(66,66),DV<66,66),  FV(66 ) . FVP(66 ) , 

1  FX(66),FXM(66),FY(66),FYM(66),PT(66),QT(66) 
C0MM0N/IHDX/NF.NFMAX,NP,NRH0,NGAM,L1,L2,L5.M1,M2,M3, 

IIST, JST, ITER, LAST, RELAX! 13), TIME. DT. XL. YL. 

2IPREF, JPREF,LSOLVE(10),LPRINT<13),LBLK(10),MODE,NTIMES(10),RHOCON 
COMMOH/CNTL/LSTOP, LSHICH 
COMMON/SORC/SMAX, SSUM 
COMMON/COEF/FLOU.DIFF.ACOF 
C0MM0N/TRAMS/CTIME,DELT,TLAST 

COMMON/OLD/  UOLOC 66 , 66 ) , VOLOC 66 ,66 ) . EOLD(66 , 66 ) , IB0DY(66 ,66 ) 
DIMENSION  U(66,66) ,V(66,66),PC(66,66),C00L(66,66) 

DIMENSION  T1(66,66),EN(66.66} 

EQUIVALENCE(F(1,1,1),U(1,I)),(F(1,1.2),V(1,1)),(F(1,1,3),PC(1,1)) 
DIMENSION  TH(66),THU(66),THDIF(66),THCV(66),THCVS(66),ST0(66,66) 
EQUIVALEHCE(X,TH),<XU,THU),(XDIF,THDIF),(XCV,THCV) 

1  ,(XCVS,THCVS).(XL,THL) 

C»KX9(KXI(Ki(XXXXXK9(Kl(XKXXi()tK»9(l(Xi(KX«XXKifX)()(XKXKX9fK)C)(KX)(KX}(XXK}CXK)(KXKi(i(9l 
DIMENSION  E(66,66),T(66,66),C0ND(66,66), 
$T0LD(66,66),RH00LD(66,66) 

EQUIVALENCE  ( F( 1 , 1 , A ) , E( 1 , 1 ) ) 

EQUIVALENCE  ( F{ 1 , 1 , 5 ) , T( 1 , 1 )  ) 

EQUIVALENCE  ( F< 1 , 1 , 10 ) ,COND( I , 1 ) ) 

C 

ENTRY  GRID 

GRID  PARAMETERS 

Hl=0.0078 

WIDTHsO.023 

DEPTH=0.5XH1 

N1X=A 

N2X»12 

N3X=A 

NAX=A 

N5X=16 

N6X:8 

N1Y=6 

N2Y=8 

N3Y=12 

NAY=12 

N5Y=18 

AL1=0.5XH1 

AL2=3.XH1 

AL3=1 .XHl 

ALA=0.5XHI 

AL5=A.5XH1 

AL6=0.5*H1 

BL1=2.XHI 

BL2=2.XH1 

BL3=1 .XHl 

BLA=2.xh1 

BL5=3.XH1 

DX1=AL1/N1X 

DX2=AL2/N2X 

DX3=AL3/N3X 

DXA=ALA/NAX 

DX5=AL5/N5X 

DX6=AL6/N6X 

DY1=BL1/N1Y 

DY2=BL2/N2Y 

DY3=flL3/N3y 

DYA=BLA/NAY 

DY5=BL5/N5Y 

XU(2)=0.0 

LFIRST=3 
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LLAST=NlX+2 
DO  10  I=LFIRST,LLAST 
XU(I)=XU(I-1)+DX1 
10  CONTINUE 

LFIRST=LLAST+1 
LLAST=LLAST+N2X 
DO  20  I=LFIRST,LLAST 
XUCI)=XU(I-11+DX2 
20  CONTINUE 

LFIRST=LLAST+1 
LLAST=LLAST+N3X 
DO  30  I=LFIRST,LLAST 
XU(I)=XUC:-1)+DX3 
30  CONTINUE 

LFIRST=LLAST+1 
LLAST=LLAST+NAX 
DO  40  I=LFIRST,LLAST 
XU(I)=XU(1-1)+DX4 
40  CONTINUE 

LFIRST=LLAST+1 
LLA5T=LLAST+N5X 
DO  50  I=LFIRST.LLAST 
XU(I)=XU(I-1)+DX5 
50  CONTINUE 

LFIRST=LLAST+1 
LLAST=LLAST+N6X 
DO  60  I=LFIRST,LLAST 
XU(I)=XU(I-1)+DX6 
60  CONTINUE 
L1=LLAST 
XLsXU(LLAST) 

YV(2)=0. 

LFIRST=3 
LLASTsNlY+2 
DO  II  I=LFIRST,LLAST 
YV(I)sYV<l-l)+DYl 
11  CONTINUE 

LFIRST=LLAST+1 
LLAST*LLAST+N2Y 
DO  21  I=LFIRST,LLAST 
YV<I)=YV(I-1)+DY2 
21  CONTINUE 

LFIRST=LLAST+1 
LLAST=tLAST+N3Y 
DO  31  I=LFIRST.LLAST 
YV<I)*YV(I-1)+DY3 
31  CONTINUE 

LFIRST=LLAST+1 
LLAST=LLAST+N4Y 
DO  41  I=LFIRST,LLAST 
YV(I)=YV(I-1)+DY4 
41  CONTINUE 

LFIRST=LLAST+1 
LLAST=LLAST+N5Y 
DO  51  1=LFIRST, LEAST 
YV(I)sYVCI-l)+0Y5 
51  CONTINUE 
M1=LLAST 
YL=YV(LLASTJ 
RETURN 
C 

ENTRY  START 

NHTRX»N1X+N2X+N3X+2 

NHTRY=HlY+N2Y+2 

NHTRX1=NHTRX+N4X-1 

NHTRY1=NHTRY+N3Y-1 

NB0STX*NlX+H2X+2 

NB0FNX*N1X+N2X+N3X+N4X+1 

NB0STY=NlY+2 

NB0FNY=N1Y+N2Y+N3Y+N4Y+1 

PRINT  «, 'INPUT  LAST  TIME  VALUE  IN  SEC. 

PRINT  K,'TIME  STEP  INTERVAL' 


AMD* 
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non 


READ  x.TLAST.DELT 
PERI=20.0 
DT=DELT 
CTIME=0. 0 
DTIME=0.0 

PRINT  *,  'INPUT  DESIRED  TIME  FOR  CONTOUR  PLOTING  IN  SEC.  ’ 
READ  x.TRIMT 
REl.IUD  (8) 

READ(8)  LTM 
DO  3000  1=1. LI 

3000  READ  (2)  ( U{ 1 , J ) , J=1 , Ml ) 

DO  3001  1=1, LI 

3001  READ  (8)  ( VC  I . J ) , J  =  1 , Ml ) 

DO  3002  1=1, LI 

3002  READ  (ECI.J),J  =  1  Ml) 

DO  3003  1=1, LI 

3003  READ  C8)  ( F( I . J , 5 ) , J=1 , “1 ) 

DO  3004  1=1, LI 

3004  READ  to-  (Pf I, J), J=1,M1) 

DO  3005  1=1, LI 

3005  READ  <8)  <  STOC I .  J ) ,  J=:l ,  Ml ) 

CP=^4178  . 

CPS=1500. 

CPHTR=439. 

RHOF=995. 

RH0S=1180. 

RHOH=8510 . 

BETA=319.25E-6 

G=9.807 

01=6. 

ALFA=0.6 
ALFAM=0.4 
DYtlV:S  =  774.2E-6 
C0HFLD=0.6027 
CONPLXsO.20 
C0NHTR=ll,7 
PI=3. 14159256 
PRsDYMVISxCP/CONFLD 
HRITEC*,*)  •PR=',PR 
RETURN 


ENTRY  SETTIM 

PRINT  «,  'CTIMEs', 'DTIME=',DTIME 
IFCDTIME.LT.lO.l)  THEN 
EIN=< .36+0.07«DTIME/10.)/WlDTH 
ELSE 

EIN=(. 430-(DTIME-10.)K. 07/10. )/MIDTH 
ENDIF 

ENGEN=EIM/CDEPTH«H1) 

P0W=ENGEN»WIDTHKDEPTH»H1 
IPRTIM=0 

WRITE  (1,«)  CTIME,E(NHTRX+1,NHTRY+6),P0W 
IFCCDTIME.Gl .9.95) .AND. (DTIME.lt. 10. 05))  IPRTIM=1 
IFCDTIME.GT.CPERI-.D)  THEN 
DTIME=0. 

IPRTIM=1 
ENDIF 

IF  (IPRTIM.EO.l)  THEN 
DSUM=0.0 
POHER»0.0 
POWER1=0.0 
DTOT=0.0 
00  504  1=2, L2 

P0WER=C0HFLD»(XU(I  +  l)-XU(I))«E<I,2)/C(YVC3)-YV(2))/2.  ) 
P0WERl=C0NFLD«(XU<I+I)-XU(I))KE(I,M2)/((YV(Ml)-YV(M2))/2. ) 
DSUM=P0WER+P0WER1+DSUM 
504  CONTINUE 
DSUM1=0.0 
POWER2=0.0 
POWER3=0.0 
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ou 


506 


999 


°°  P§WER2=C0HFLD5«(YV(J+1)-YV{J))»E(2,J)/((XUC3)-XU(2))/2.) 

POWERS =C0HFL D* < YV( J+1 )-YV( J ) )XE(L2»J)/(CXU(L1 )-XU( L2) )/2 . ) 
D5UM1=DSUM1+P0HER2+P0HER3 
CONTINUE  - 

0UT=' ,DSUM+D5UM1 


•ENERGY 

•ENERGY 


WRITE!*.*) 

WRITE!*.*)  _  - - 

WRITE!*.*)  •CURRENT  TIM=^.CTIME 
ENDIF 

FARK=TRINT-CTIME 
IF  !TRINT.EQ.CTIME)  THEN 
DO  999  1  =  1. A 
LPRINT!I)=.TRUE. 

CALL  PRINT 
LPRINT!I)=. FALSE. 

CONTINUE 
NTIME=1 

WRITE!6.«)  •START  WRITING  ON  UNIT  9' 
WRITE  !*.*)  •ONCE  YAZDIM* 

REWIND  !3) 

WRITE  !3)  NTIME.Ml.Ll 
WRITE  (3)  !YV!J).J=2.M1) 

WRITE  !3)  !XU!I).I=2.L1) 

WRITE  !3)  (YCJ). J=1.M1) 

WRITE  (3)  !X!I).I=1.L1) 

WRITE  (3)  CTIME 
DO  300  J=I.M1 

300  WRITE  (3)  (V(I,J),I=1,L1) 

DO  301  3=I>MI 

301  WRITE  (3)  <U!I,J),I=1.L1) 

DO  302  J=1.MI 

302  WRITE  <3)  <E<I. J),I=1.L1) 

DO  303  J=1.M1 

303  WRITE  (3)  (P(I, J),I=1,L1) 

DO  305  Jsl.Ml 

305  WRITE  <3)  (STO(I,J),lsl.Ll) 

DO  304  Jsl.Ml 

304  WRITE  C3)  <F(I, J,3),1*1,L1) 

DO  306  3=1. Ml 

306  WRITE  (3)  (IB0DY!I.J).I»1.L1) 
WRITE(6.«)  'CLOSE  UNIT  3* 

CLOSE  (3) 

ENDIF 

C  STORE  OLD  VALUES 
DO  081  1=1. LI 
DO  881  J=1>M1 
UOLD!I.J)=U(I.J) 

VOLOCI.J)=V(I,J) 

EOLD!I.J)=E<I,J) 

881  CONTINUE 

LSTOP=. FALSE. 

LSWICH=. FALSE. 

ITER=0.0 

DTIME=DTIME+DELT 

RETURN 


C 

C 


501 

500 


502 


ENTRY  DENSE 
RETURN 

ENTRY  BOUND 
RETURN 

ENTRY  OUTPUT 
IFIITER.NE.O)  GOTO  500 

FORMAT! '^^^^ITER'.BX, 'SMAX'.BX,'  SSUM'.TX. 'VEL! FRONT) 
16X. •T(HEATER)',6X, 'CURRENT  TIME’) 

CONTINUE 

WRITE! *. 502)  ITER, SMAX. SSUM, V(NBOFMX+l , NHTRY+6 ) 

1 , E! NHTRX+1 , NHTRY+6 ) , CTIME 
FORMAT  (I6.1P5E12.3) 
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400 

401 

402 

403 

404 

405 


RETURN 


ENTRY  GAMSOR 
IF  (NF.EQ.3)  RETURN 
IF  (NF.EO.l .0R.NF.E0.2)  THEN 
DO  400 

DO  400  I=1>L1 
RH0(I,J)=RH0F 
GAM(I.J)=DYNVIS 
CONCI, J)=0.0 
IBODYCI, J)=l 
CONTINUE 

DO  401  I=NBOSTX,NBOFNX 
DO  401  J=NBOSTY,NBOFNY 
GAM(I,J)=1 .E+20 
eONTINUE 
ENDIF 

IF(NF.EQ.2)  THEN 
DO  402  1=2, L2 

TEM=(FY(J)»E(I, J))  +  (FYM(J)«E(I, J-D) 
C0N1=RH0F«GJ«BETA»(TEM-TSTART) 

C0N( I , J ) = ALFAKCONl+AL  FAM»STO( I , J ) 
STO(l,J)=COHl 
CONTINUE 
ENDIF 

if{nf.eo.4)  then 
DO  403  J=1,M1 
DO  403  1=1, LI 
RHOCl, J)=RH0F»CP 
GAM(I,J)=CONFLD 
CON(I,J)=0.0 

continue 

DO  404  I*NB0STX,HB0FNX 
DO  404  J=NBOSTY,NBOFNY 
RH0(1,J)=RH0S«CPS 
GAM(I,J)=CONPLX 

CONTINUE 

00  405  I=NHTRX,NHTRX1 
DO  405  J=NHTRY,NHTRY1 
RHO(I,J)=RHOH*CPHTR 
GAM<I,J)=COHHTR 
CON(I,J)=ENGEN 


CONTINUE 

ENDIF 

RETURN 

END 
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APPENDIX  B.  STEADY  STATE  AND  TRANSIENT  POWER  PROGRAMS 

FOR  HP  85 


10 

1  !  I 

1  1  ,  1  t  1  1  1  1  t  1  i  1  1  i  1  1  M  1  1  1  1  1 

20 

1 

POUEP  PROGRPH  FOR 

30 

1 

HP-59501B  PROGRPMMER  ! 

40 

1 

TO  SUPPLY  STEPDV  '..'OLT ' 

50 

i 

TO  THE  HERTER  PT  THE  ' 

€’  0 

1 

294 

BOTTOM  OF  CEHT.  COL.  ! 

70 

1  1  ! 

1  1  t  C  t  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

I 


so  OISP  "INPUT  MfiX.  PND  PUE .  'JO 
LTPGE" 

SO  INPUT  ill.Ul 
100  Rl=Ul/lll»10OO+20O0 
110  OUTPUT  70iS  USING  “#,4D"  ;  R1 
120  END 


10  !  I  M  <!!  M  I  I  I  I  I  I  I 

!  . 

20  1.  TRANSIENT  POWER  PROGRRPM  ' 
! 

30  !  f  !  !  !  >  M  !  ■  ■  ■ < ! <  > I ■ ) >  I  I  I  I  I  I  I 

I  I  . 


40  < 

50  DIM  E<:  I??) 

SO  OISP  "INPUT  I1PX.  VOL" 
70  OISP  "PND  MIN,  VOL" 
80  INPUT  R1,R2 
90  DISP  "INPUT  PERIOD" 
100  INPUT  P 
110  0*0^ 

120  Z=iP/'C 
130  N=INTT2) 

140  K»Nx2 


150  Vl<K)*INT<Rl/25»100)»lO+200e 
160  Vl<e)siNT<R2x25:kie8>tl0-*-2O00 
170  E<K)»R1 
180  E<0>aR2 


190  L»K-l 

200  FOR  1=1  TO  L 

210  E<I)»(Rl-R2>/K+E<I-l) 

220  E<N-n=E<I) 

230  Vl<I)=INT<E<I>/25*100>*10+20 


240  Vl<N-i:'»Vl<I)  • 

250  NEXT  I 
260  M*N-l 
270  BEEP 

280  FOR  1=0  TO  H 

290  OUTPUT  706  USING  "§.40"  ;  VI 

<o 

300  WRIT  56 
310  NEXT  I 
320  GOTO  270 
330  END 
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APPENDIX  C.  TEMPERATURE  AND  POWER  ACQUISITION  ROGRAM 


FOR  THE  HEATER  AT  THE  BOTTOM  OF  THE  CENTERAL  COLUMN 

10  I  I  I  I  I  M  I  I  I  I  I  M  M  I  I  I  I  I  I  I  I  I  M  I  I  I  !  I  I  j  j  !  I  I  M  M  *  I  I  !'>  I  M  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

20  Ml  TRANSIENT  ! i ! 

30  Ml  temperature  and  power  acquisition  program  for  the  heater  I  I  I 

40  ! ! I  AT  THE  BOTTOM  OF  CENTERAL  COLUMN  i i i 

50  I  I  I  I  I  1  I  I  1  I  I  I  I  I  I  I  I  M  I  I  I  I  I  I  I  1 ! I M I !! 11  I  I M ! j  I  M  1  I  1  I  I  I  I  I  I  I  I  I  I  I  I  1  M  I  I  I 

60  OPTION  BASE  I 

70  REAL  TuU  350  )  ,Yig(  180  )  ,Cas<350  )  ,Tei>i<  180  )  .UoK  350  )  .Ba t h(  3  )  .Bat  (  3  ) 

80  CREATE  BOAT  " Y I6I T :  ,700 , 1 ■ . 1 80 
90  ASSIGN  @Path4  TO  " Y IG I T :  ,700 , 1 " 

100  CREATE  BOAT  "CANSU:  .700 . 1 ' .350 
110  ASSIGN  3Path5  TO  "CANSU :  ,700 . 1 " 

1 20  Res i 5 t = 1 . 886 

130  R2“Re5ist»Resist 
140  BEEP 

150  FOR  I«1  TO  350 

160  OUTPUT  709; "CONFMEAS  DCU  ,508-509  .USE  0" 

170  ENTER  709;U,U1 
180  Uol(I)=U 
190  TuItI)»U1 
200  NEXT  I 
210  FOR  1-1  TO  180 

220  OUTPUT  709; "CONFMEAS  DCU  .209 .USE  0" 

230  ENTER  709;T 

240  WAIT  .09 

250  Ten(l)-T 

260  NEXT  I 

270  FOR  1-1  TO  350 

280  IF  I>180  THEN  GOTO  30® 

290  Yig(  I  0006797+< 25825. 1328»Tei<i(  I  >)-< 607789. 2467»(  Ten(  I  )»Ten(  I  ) )  >-(2195203 
4.33G4»<Ten( I  )"3  )  >+(8370810996. 1874- (Tent  I )"4) >+.039 

300  Ca8( I >-(Uol( I >-Tul( I >  >/2.0l»(Tul<I )-<Uol( I  )-Tul( I >  >/2 .01 •(  4700/ 1 47. >  > 

310  NEXT  I 

320  OUTPUT  709; “CONFMEAS  DCU  .31 7-31 9  .USE  0* 

330  FOR  I-I  TO  3 
340  ENTER  709;Bat( I > 

350  NEXT  I 
360  FOR  1-0  TO  3 

370  Bath( I >-.000B7g7+( 25825. 1328-Bat < I > >-(607789. 2467* (  Bat ( I >*Bat(  I >  >  >-(219520 

34.3364*(Bat(  I  >  “3  >>+(  837081 0996 . 1874*(Bat(  I>'‘4  )  >+.036 

380  NEXT  I 

390  A$-‘TOP(C>' 

400  B$-"MIO(C)' 

410  C$-"B0T(C>" 

420  PRINT  USING  •2X  .AAAAAAAAAAAAAAAAAAAAAAA .3X .00.00" ; AS ,Bath(  11 
430  PRINT  USING  ■2X  .AAAAAAAAAAAAAAAAAAAAAAA  .3X  .00 .00" ;BS  .Bat h( 2 > 

440  PRINT  USING  ■2X  .AAAAAAAAAAAAAAAAAAAAAAA  .3X .00. DO" ; CS  .Bat h( 3 > 

450  PRINTER  IS  CRT 
460  OUTPUT  8PathS;Ca8( •  ) 

470  OUTPUT  @Path4;Yig(* > 

480  END 
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APPENDIX  D. 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
I  10 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 


TEMPERATURE  ACQUISITION  PROGRAM 

I  I  I  I  I  I  I  I  I  M  M  I  I  I  I  I  !  !  1  !  !  !  M  !  !  !  I  M  I  I  I  !  !  !  M  !  !  I  !  M  I  !  t  I  !  I  '  !  I  I  I  I  M  M  !  !  I  I  !  !  I 

temperature  ACQUISITION  PROGRAM  " 

FOR  15  HEATERS  IN  (STEADY  STATE)  N 

I  ,  I  M  I  I  I  I  I  I  I  J  !  I  !  !  !  !  I  !  t  !  M  M  !  !  !  I  I  M  !  !  M  !  !  I  I  I  I  I  I  I  I  !  I  I  M  I  I  I  M  I  I  I  M  I  I  r  M  I 

OPTION  BASE  1 
REAL  Uolts< 18  ) 
real  Tenp( 18) 

PRINT  DATES(TIMEOATE) 

PRINT  “INPUT  AUERA6E  POWER' 

INPUT  Power 
PRINTER  IS  701 
PRINT  OATES( TIMEOATE ) 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT  "  “ 

PRINT  "  ” 

PRINT  '  “ 

I  I  I  I  I  I  (I  I  I  I f I  I M I n  M  M ' H 11  SURFACE  f I -5 
OUTPUT  709i “CONFMEAS  DCU ,  115-1 19, USE  0" 
for  1-1  to  5  ^ 

Tei>ip<  I  )«  .000G797+<  2S82S.  t328»Uolt5<  I  )  )-<  607789. 2467»(  UoU  si  I  >»UoHs(  I  ) ) )-( 


Power,  “WATT  INPUT  POWER' 
15  HEATER  ELEMENTS  ' 

3  MM  SHROUD  SPACING 


21 952034. 3364. < Uol ts< I )'3))+< 8370010996. 1874. (Uolts(  I  )'4))+. 039 


250 

260 

270 

280 

290 

300 


next  I  .  ,,, 

I  I  M  I  I  I  I  I  I  I  M  M  !  M  M  M  1  M  M  SURFACE  t2 1  -#30 
OUTPUT  709i 'CONFMEAS  OCV ,200-209 .USE  0' 

FOR  I-G  TO  IS 

TeIpaI-!00067g7+<2S82S.l328.Volts(I >>-<607789. 2467. (Volt5( I LUoltsI I > > )-< 


2ig52034.3364.(VoIt9<I>'3>>+<83708l099G.I874.(Volt5(I >'4>)+.039 


310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 


NEXT  I 

1  !  !  PRINT  OUTPUT  TEMPERATURES 
PRINT  •  • 

FOR  1-1  TO  IS 
PRINT  "HEATER  •'  ll 
NEXT  I 
PRINT  ■  ' 

PRINT  '  ' 

PRINT  ■  ' 

PRINT  ■ 

PRINT  ■  ' 

OUTPUT  709; ’CONFMEAS 
FOR  1-16  TO  18 

Tenp”n-!0006797+(  25825. 1328.Volts(  I  >  >-(  607789. 24E7.(Uolts(  I  LUoltsI  !>>>-( 


TEMPERATURE  IN  DE6.  C  -•;Tenip(I> 

BATH  TEMPERATURES  (TOP  TO  BOTTOM)" 
DCV,3t7-3l9,USE  0" 


Hail  - - - -  -- 

2l952034.3364.<Uolt9(  I )*3 >)♦< 8370810996. 1  874. ( Uol 1 5( I  )'4 > >+.036 


460  PRINT  "TEMP.  DEG.  C  ■;Tei>ip(I> 

470  NEXT  I 

480  PRINTER  IS  CRT 

490  ENO 
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APPENDIX  E.  POWER  INPUT  PROGRAM 

40  OUTPUT  709) "RST” 

50  OPTION  BASE  t 

B0  REAL  Uolts( 18 ) 

70  REAL  Pou( 15) 

80  REAL  Ptot  ,Pavg  .Perdi F 

'90  OUTPUT  709; “CONFMEAS  DCU.3IG,USE  0" 

100  ENTER  709;Uolt9( 16  ) 

1  10  OUTPUT  709) "RST" 

120  OUTPUT  709; "CONFMEAS  DCU  .4 1 0-4 1 7  .USE  0" 

130  FOR  I-l  TO  8 

140  ENTER  709;Uolts( I ) 

150  NEXT  I 

160  OUTPUT  709; "RST" 

170  OUTPUT  709; "CONFMEAS  OCU.401,USE  0" 

180  ENTER  709:Uoit5( 18 ) 

190  OUTPUT  709;"RST" 

200  OUTPUT  709) "CONFMEAS  DCU  .503-509 .USE  0" 

210  FOR  1=9  TO  15 

220  ENTER  709;UoIts( I ) 

230  NEXT  I 

240  OUTPUT  709; "RST" 

250  OUTPUT  709) "CONFMEAS  DCU  ,500  .USE  0" 

260  ENTER  709)Uoas<  17) 

270  OUTPUT  709) "RST" 

280  PRINT  DATES! TIMEOATE ) 

290  PRINT  "  " 

300  PRINT  •  ■ 

310  PRINT  "INPUT  UOLTAGE  D.C.  UOLTS  PAN  f;Volt9(18) 

320  PRINT  "INPUT  UOLTAGE  D.C.  UOLTS  PAN  2")Uolts(17) 

330  PRINT  •  " 

340  PRINT  •  ' 

350  Re9i9t‘‘2.0 

360  R2“Resi9t»Re9i9l 

370  Pow(3)-<Uolt9( 18)-Uolt9(3))«(Uolt9(3)»<Re9i9t+.55)-Volt9( 18)».55)/R2 

380  Pou(8  >-<UoU9(  18  )-Uolt9<8  )  )•<  VoltalB  )*<Re9i9t  +  .55  )-Uolt9(  18  )».55  )/R2 

390  Pou( 10  )-(UoIt9( l7)-UoIt9< 10) »*(Volt9< 10 )*< Re9 19 1  + . 3 )-Uol 1 9< I8)».3  )/R2 

400  Re9i9t-1.972 

410  R2=Resi9t»Real9t 

420  Re9 i9t=1 . 886 

430  R2=Re9i9t*Re9l9t 

440  Poui(  rS)-<UoIt9<  I4)-Uolt9(  l5))/2.0»»<Volts(  l5)-((Volt9(  U)-Volt9(  I5>>/2.0I» 
4700/147. ) ) 

450  Re9i9t”l .862 

460  R2*Re5i9t»Re9i9t 

470  Pou< 12  )-(Uolt9( 17)-Uolt9( 12  ))•( Volts! 1 2  >•( Res i s t  + . 5 )-Uol t s< 17 )».5 )/R2 
480  Pou( l5)-<Uolts< I7)-Uolts( 13  ))•< Volts! J3)»( Resist +0. >-Uo]ts! I7)»0.  )/R2 
490  Resist"! .852 

500  R2”Re9ist»Re9l9t 

510  Pou!9)-!Volt9!17)-VoUs!9>)»!Volts(9>»!Re9l9t  +  .l  )-Volts!  17)».1  )/R2 

520  Resist"! .8182 

530  R2"Resi9t»Re9i5t 

540  Pou!6  )"!Volt9! 18  )-Volt9!6  )  )•! Vo  1 1 s! 6 )•! Resist  + . 7 )-Vol t s!  17 )».7  )/R2 

550  Resist"! .8 

560  R2"Re9ist*Re9i9t 

570  Pow!4 )"!  Volts! 18 )-Volts!4 ) )*!Volts!4 )»!Re9lst  +  .2  l-Volts! 18  )».2  )/R2 
580  Pow!5)"!Volt5! 1 8  )-Vol t s! 5  )  )*! Vol ts! 5  )•! Resist*.? )-Vol t s! 18 )«.7 )/R2 
590  Pow! 14  )"!Pom! 15  )+Pou! 13) )/2. 

R00  Poui!  7  )-!Pou!  6  )*Pou!  8  )  )/2 
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GIB  Re5i5t=I. BEE7 
G20  R2=Re5lst ‘fiesj  st 

G30  Pou)(  1  )  =  <  Uol  ts(  1  8  )-Uol  ts(  1  )  )•(  VoltsC  I  )•  (  Res i  s t  +  .  4  )-iJd1  t s <  1  8  )•  .4  )/R 

640  Pou(  2  )-('Jolt5(  I  8  >-iJoUs(2  )>»(yolts(2)o(Resist+.4  )-Uol  ts(  1  8  )•  .4  )/R 

650  Poui(  n  )•(  yol  ts<  17  )-Uoll5(  n  >  )•(  Uolts(  l  l  )*(Resi5t+.2  )-VoUs(  17  )•  .2  )/R2 
6G0  FOR  1=1  TO  15 

670  PRINT  "BLOCK  #  "  i  I  . "  POUER  -  "jPoud)."  UftTTS" 

680  WAIT  .5 

630  NEXT  I 

700  Pto4=0 

710  FOR  1=1  TO  15 

720  IF  1=14  THEN  740 

730  Ptot=Ptot+Pou( I  ) 

740  NEXT  I 

750  Pavg=Ptot/l4 
760  PRINT  ■  " 

770  PRINT  "AOERAGE  POUER  -  "iPavg."  WftTTS" 

780  PRINT  "  ■ 

790  FOR  1=1  TO  15 

800  Perdi f-ABS(  (  Pavg-Pou( I ) )/Pavg )«100 

810  PRINT  "HEATER  fil."  PERCENT  OEOIATION  "iPerdif 

820  WAIT  .5 

830  NEXT  I 

840  END 
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APPENDIX  F.  TEMPERATURE  ACQUISITION  AND  REDUCTION 

PROGRAM  (STEADY  STATE) 

I  I  I  I  I  I  I  I  I  I  L  I  I  M  M  I  I  I  I  I  I  I  I  I  I  t  I  I  I  M  I  I  I  I  f  I  I  I  I  I  I  f  M  I  M  I  I  I  f  I  I  I  M  I  (  I  I  M  I  I  I  M  I 

20  ii  TEMPERATURE  ACQUISITION  AND  REDUCTION  PROGRAM  m 

30  !l  STEADY  STATE  (CENTER  COL)  I' 

40  ! I  PROGRAM  PERFORMS  THE  REQUIRED  FUNCTIONS  NEEDED  TO  OBTAIN  DATA  ' ' 

50  II  FROM  THE  ACQUISITION  SYSTEM  AND  CONVERT  THE  THERMOCOUPLE  OUT-  II 

B0  II  PUT  VOLTAGES  INTO  TEMPERATURE  VALUES.  ii 

70  M  M  I  I  I  I  I  I  I  I  I  I  I  M  I  I  M  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  M  I  I  I  I  I  I  I  I  I  I  M  I  I  I  I  I  I  I  I 

80  OPTION  BASE  1 

30  REAL  Vo  It s( 18) 

100  REAL  T< l8),Xtol( l5).Nu( l5>.Ra( I5),0elt(l5) 

1 10  REAL  Lnu( IS ) .LgrC 15 ) 

120  REAL  Power  ,Tanb  ,G  ,Kpg  .Oelx  , As  .Per  .Ta  ,Oconv  ,Lba'" 

130  REAL  Of  lux  ,T f 1  In  .Beta .Spvol  .Oynvisc .Ocond 
140  REAL  K 1  nvi  sc  .K  f  .Pr  .The  t a  .H.x  .Grn  .Hux 
150  PRINTER  IS  CRT 
160  PRINT 

170  PRINT  "ENTER  AVERAGE  POWER  -WATTS" 

180  INPUT  Power 

190  PRINT 

200  PRINT  "ENTER  SHROUD  SPACING  IN  nn" 

210  INPUT  Spaces 

220  iPRINT  •  " 

230  I  M  I  I  I  I  M  '  I  I  I  M  I  I!  I  !  M  I  !  M  II  SURFACE  #16-20 
240  OUTPUT  709! "CONFMEAS  DCV ,  I  15-119, USE  0* 

250  FOR  I-l  TO  5 
260  ENTER  709!Volts< I > 

270  T(  I )-. 0006797+( 25825. 1328»VoUs( I ))-( 607789. 2467* (Volt s( I >»Volt5( I  ) )  >-(219 
52034. 3364* ( Volt s( I T'S > )+( 83708 10996 . 1874»( Vol t s( I >"4 ) )+ .039 
280  XtoK I  )-(31 .75-6.g85-( I-l 6  )• 1 .27  )/31 .75 
290  NEXT  I 

300  I  I  I  I  I  I  I !  I  I  I  I ! (  I  M  M  I  1 M  M  M  SURFACE  i21-#30 
310  OUTPUT  709!"CONFMEAS  DCV  ,200-207  .USE  0* 

320  FOR  1-6  TO  13 

350  IF  1-14  THEN  370 

340  ENTER  70giVolts(I) 

550  T( I  )-.000G7g7+( 25825. 1328»Vol 1 8(1 )  )-(607789.2467«( Volts( I )«Volts( I) ) >-(219 
52034. 33G4»( Volt s( I  )  "3  ))+( 83708 1 0996 . 1874»{ Vol ts( I  )'4))+.059 
360  XtoK I  )-(3l ,75-G.g85-( I- 1 6 )* 1 .27 )/31 .75 
370  NEXT  I 

380  OUTPUT  709; "CONFMEAS  DCV  ,209  .USE  0" 

390  ENTER  703iVolts( IS) 

400  T(  15  )-.000G797+(  25825. 1328*VoUs(  15  )  >-(  607789. 2467»(  Volt  s(  1 5  )»Vol  t  s(  1 5  ) )  )- 
(21952034.3364M Volts( 15)“3 ) )  +  ( 8370810936. 1 874»( Volt s( 15 )"4  )  )+.090 
410  OUTPUT  709!"CONFMEAS  DCV  ,3 1 7-31 9  .USE  0" 

420  FOR  I-I6  TO  18 

430  ENTER  709(Volts( I  ) 

440  T( I  )-.000B7g7+( 25825. 1328»Volts( I > )-( 607789. 2467»( Vol ts( I )»Volts( I ) ) >-(219 
52034. 3364»( Volt s< I  )'3  ))  +  ( 8370810996 . 1874»( Volts! I  )'4  )  >+.036 
450  PRINT  "TEMP.  DEG.  C  "iTd) 

460  NEXT  I 

470  PRINT  "INPUT  AMBIENT  TEMPERATURE" 

480  Tanb-< T( 16  )+T(l7)+T( 18  )  )/3 
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490  PRINTER  IS  701 
500  PRINT  “  ■ 

510  PRINT  DATE$(TIMEDATE) 

520  PRINT  "DATA  REDUCTION  OUTPUT" 

530  PRINT  "  " 

540  PRINT  "POWER  LEUEL  (WATTS)  :"  .Power 

550  PRINT 

560  PRINT  "SHROUD  POSITION  ( nn  )  Space* 

570  PRINT 

560  PRINT  "AMBIENT  TEMPERATURE  (C)  :”.Tanb 

590  PRINT 

600  PRINT 

610  I 

620  I 

630  6-9.81 

640  Kpg= .1421 

650  Deix-. 006751 

660  As-. 0239* .0070 

670  Per=2* .0239+2* .0078 

680  Lbar-As/Per 

690  A$="0cond" 

700  B$=")(/L" 

710  D$="L06(Nu.<  )" 

720  ES="L06(6r.y*  )" 

730  FS--OELT" 

740  65- "H>," 

750  H$-"Nu-;" 

760  IS="6r.r»" 

770  J$-"Tenp" 

780  PRINTER  IS  CRT 

790  I  M  I  '  I  I  I  M  I  I  I  I  I ( ! M I  I  1  I  '  I  I M ! > I M  M ! M  I  I  !  !  I  I  I  !  M 
P00  !  ENTER  CONOUCTION  LOSSES  1 

610  !  I  M  I  1  1  I  1  I  H  M  I  I  I  >  '  '  !  1  M  M  I  ■  H  I  M  M  !  M  M  !  I  *  M  I  I  I 

620  PRINT  "INPUT  CONDUCTION  LOSSES  AS  COMPUTED  FROM  ELI.PACK" 

630  INPUT  Ocond 

640  M  !  I  I  I  I  I  !  I  I  M  11  !  n  I  I  I  M  i  I  I  I  !  I  I  I  I  !  M  I  I  I  I  I  M  I  !  I  !  !  ! 

850  PRINTER  IS  701 

860  PRINT  "  CENTER  COLUMN" 

870  PRINT  "  " 

880  FOR  I-l  TO  15 

890  PRINT  "  " 

900  PRINT  "  SURFACE  t" ,I 

910  PRINT  ■  " 

920  PRINT  USIN6  " 2X  , AAAAA  ,5X  .AAA  .3X  .AAAAAAAA  .2X  .AAAAAAAAA  ,3X  ,AAAA ,4X  . AA  .4X  .AAA 
.3X  .AAAA.2X  .AAAA" ;  A$  .8$  .DS  .ES  .F$  .6$  .H$  .1$  ,J$ 

930  Oconv-Pouei — Qcond 

940  Of  lu.<-0conv/As 

950  Delt( I  )»T(I  )-Tanb 

960  Tfiln-273. 15+(T(I  )+Tanb)/2 

970  Beta-(  -7 . 944858E-8* ( T f i ln"2 ) >+( 5.7356479E-5»Tf i In  >-.0097610563 
980  5pvol-(4.69699E-9»<Tf iln'2  )  )  +  < -2 . 53745E-6»T f i In  >+.001341903 

990  Dynvi5C-( 3. 234851 1E-7«(Tfiln"2 )  )+(-.0002l4744B7*Tf iln  >+.036166792 
1000  Kinvisc“5pvol*Dynvisc 
1010  Kf-(1 .41 8181 818E-3*Tf iln >+.1866 

1020  Pr-(4.6570BE-3*(TfiIn'2))+(-2.922094»Tfiln>+463.3319  • 

1030  Thetal-Pouerl /( .0220»Kf > 

1040  Tenp-T(I> 

1050  Theta-Qf luy*Lbar/Kf 
1060  Hx-Of 1UA/Delt( I > 

1070  6rn-( G*Bete*Qf lux»< Lbar  >"4  )/(Kf»Kinvisc“2  > 

1080  Nux-Hx*Lbar/Kf 
1090  Nu( I >-Nux 
1  1 00  Ra(  I  >-Pr*Gri^ 
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1110  Lgr( I  )=LGT( Grn  ) 

1120  Lnu(  I  )“LGT(Nuy,  ) 

1  I  30  LQ  =  Ugr( I  ) 

1 1 40  Lnn»Lnu( I > 

1150  PRINT  USING  "  I  X  .0 . DODD  .2X  .0 . ODD  .2X  .0 .  ODD  ,5X  .D .  ODD  .4X  .DO .  ODD  ,2X  .DDDD .  D  .  1  X  .0 
O.ODD  .DDDDD.D  ,1  X  ,DD  .DD"  iQcond  ,Xtol(  I  )  .Lnn  .Lg  ,DeJt<  1  >  ,Nu)>  ,6rn  .Temp 

1160  NEXT  I 

1170  CREATE  BOAT  "NU:  ,700  ,0  ,2'' ,  15 
1180  ASSIGN  @Path4  TO  "NU: .700  ,0 .2 " 

1190  OUTPUT  ePathAjNufO 

1200  CREATE  BOAT  "RA :  .700 , 0  .2 “  ,  1 S 

1210  ASSIGN  @Path5  TO  "RA:  .700  ,0  .2" 

1220  OUTF JT  ePath5iRa( •  ) 

1230  PRINTER  IS  CRT 

1240  CREATE  BOAT  "OTEMPCWSMA ;  ,700 . 1 "  .  1 5 
1250  ASSIGN  ePath7  TO  "DTEMP2U3MN : ,700 . 1 " 

12B0  OUTPUT  @Fath7iDelt< •  ) 

1270  END 
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